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Abstract. We study triangulation schemes for the joint kernel of a diffusion process with 
uniformly continuous coefficients and an adapted, non-resonant Abelian process. The proto- 
typical example of Abelian process to which our methods apply is given by stochastic integrals 
with uniformly continuous coefficients. The range of applicability includes also a broader class 
of processes of practical relevance, such as the sup process and certain discrete time summa- 
tions we discuss. 

We discretize the space coordinate in uniform steps and assume that time is either con- 
tinuous or finely discretized as in a fully explicit Euler method and the Courant condition is 
satisfied. We show that the Fourier transform of the joint kernel of a diffusion and a stochastic 
integral converges in a uniform graph norm associated to the Markov generator. Convergence 
also implies smoothness properties for the Fourier transform of the joint kernel. Stochastic 
integrals are straightforward to define for finite triangulations and the convergence result gives 
a new and entirely constructive way of defining stochastic integrals in the continuum. The 
method relies on a reinterpretation and extension of the classic theorems by Feynman-Kac, 
Girsanov, Ito and Cameron-Martin, which are also reobtained. 

We make use of a path-wise analysis without relying on a probabilistic interpretation. The 
Fourier representation is needed to regularize the hypo-elliptic character of the joint process 
of a diffusion and an adapted stochastic integral. The argument extends as long as the Fourier 
analysis framework can be generalized. This condition leads to the notion of non-resonant 
Abelian process. 
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1. Introduction 
Consider a diffusion defined by the Markov generator 

(1.1) 

on the bounded interval X = [—L^.L^] C M where < L^^ < oo. For simplicity, we assume 
periodic boundary conditions and identify the two boundary points ±Lx with each other. The 
coefficients cr(x)'^ and are assumed to be at least uniformly continuous and our results will 
depend on the degree of smoothness. We also assume that <j{x) > Sq for some constant Eq > 0. 
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In this paper, we consider path dependent processes such as stochastic integrals of the form 



where a(x) and b{x) are functions on A — [^L^, L^]. The case a{x) — was considered in 
(Girsanov 1960), (Cameron and Martin 1949), (Feynman 1948) and (Kac 1948). The case with 
a{x) 7^ was first tackled in (Ito 1949). The proofs were streamlined by using martingale theory 
in (Doob 1953). Other references for stochastic integrals arc (Kunita and Watanabe 1967), 
(McKean 1969) and (Meyer 1974). Most of these results assume smooth coefficients. The more 
general case of uniformly continuous coefficients was considered in (Stroock and Varadhan 1969). 

We are interested in convergence properties of triangulation schemes for the joint distribu- 
tion between the underlying process Xt and the stochastic integral yt- If the diffusion process is 
approximated by a sequence of Markov chains then existence and uniqueness of stochastic inte- 
grals is not problematic. The problem is to determine conditions under which joint distributions 
between stochastic integrals and the underlying process converge in some meaningful norm in 
the limit as the triangulation becomes finer and finer, approaching the continuum limit. We 
establish convergence assuming uniform continuity in the coefficients of the diffusion and the 
stochastic integral and find convergence rates. Convergence takes place in the graph-norm of 
the Fourier transformed Markov generator. In the continuum limit, the Fourier transform of the 
joint distribution is entire analytic in the conjugate variable of the stochastic integral and is in 
the domain of all the powers of the Fourier transformed Markov generator. Consideration of the 
Fourier transform is essential in the proof as the joint distribution itself can be singular even in 
trivial special cases such as the one where a{x) = b{x) = 0. The Fourier analysis also leads to 
a new derivation of the celebrated formulas for characteristic functions of stochastic integrals in 
(Girsanov 1960), (Cameron and Martin 1949), (Feynman 1948) and (Ito 1949). Furthermore, 
analyticity implies convergence of moment formulas. 

These results are classic and well established in the continuum limit. One of the new facets 
of our derivation is that it is entirely constructive and thus instructive from a computational 
viewpoint. We make no use of compactness arguments and non-constructive measure theory 
methods. The derivation is based on pathwise analysis and a renormalization group transforma- 
tion. The construction applies to discretized operators and does not rely on the existence of a 
probabilistic interpretation. The proof of convergence is actually carried out for analytic exten- 
sions of the Fourier transform of the joint distribution, so we consider complex valued weights 
for the paths. The argument would also extend to different type of diffusions with complex 
coefficients as they occur for instance in quantum mechanics. Furthermore results are extended 
to discrete time approximations with fully explicit Euler schemes, in which case we also derive 
convergence bounds in the graph- norm for Fourier transformed joint kernels. 

In (Albanese 2007c) we discuss the case of the kernel of diffusion processes, a more elementary 
situation. This paper is strictly more general and it is self-contained with no dependencies on 
our previous work as all derivations need to be adapted in detail. The key technique of pathwise 
analysis and renormalization group resummations over decorated paths are however quite similar. 

Our operator algebraic approach to the problem also points out that stochastic integrals are 
not the most general class of path dependent processes to which these methods apply. This is 
based on the recognition that Fourier transforms play a key role as they allow one to block- 
diagonalize the joint kernel. But block-diagonalization can be achieved under more general 
conditions. We consider a generic situation of a diffusion and an adapted process. We associate 
an operator algebra to a finite triangulation of the joint process. If this operator algebra is 
Abelian, i.e. commutative, and if it satisfies a certain non-resonance condition, then the joint 
Markov generator can also be block-diagonalized and methods can be extended. 

Non-resonant Abelian processes are a broad class of hypo-elliptic multidimensional diffusion 
processes which are computable by robust methodologies based on block-diagonalizations and 
fully expUcit Euler schemes. Several instances of these processes found engineering applica- 
tions already, see (Albanese and Vidler 2007), (Albanese and Trovato 2005), (Albanese 2007a), 



(1.2) 




STOCHASTIC INTEGRALS AND ABELIAN PROCESSES 



3 



(Albanese et al. 2006), (Albanese and Osseiran 2007). As an example, we discuss here the case 
of the sup processes of a one-dimensional diffusion, i.e. of 

(1.3) yt = sup Xs- 

6[0,t] 

We also discuss discrete time summations of the form 

n 

where T, = iAT. 

The kernel convergence in graph norm we establish in this paper has noteworthy implications 
from the computer science viewpoint. As we discuss in (Albanese 20076), they imply that 
explicit schemes give a robust valuation methodology for joint kernels of processes specified 
semi-parametrically. The Courant-Friedrichs-Lewi stability condition (Courant et al. 1928) for 
explicit methods is not much of an impediment in the presence of sufficient system memory 
as linear fast exponentiation can be used to accelerate the scheme. The ability of evaluating 
kernels robustly is the main reason why analytically solvable models are interesting and broadly 
used. But not all models solvable in closed form are practically computable in the sense that they 
allow for robust kernel valuations within the limits of double precision floating point arithmetics. 
Higher order hypergeometric functions are often difficult to handle and may require multiple 
precision. Linear fast exponentiation is a technique based on accelerated fully explicit schemes 
that empirically is observed to work well even in single precision and for very general model 
specifications. Empirically, we find that it performs better in single precision on some solvable 
models than the use of closed form solutions in double precision. This paper was prompted by 
the observation of this phenomenon and the desire to understand and explain the underlying 
smoothing mechanisms which is of technological significance. 

The paper is organized as follows: In Section |2] we introduce notations and state our results 
regarding stochastic integrals. Proofs in the case of stochastic integrals are in Section [3) where 
we consider the case of continuous time and in Section [4] where we discuss convergence for fully 
explicit Euler schemes. The Ito representation is in Section 5. Extensions to the sup and other 
Abelian processes are given in Section 6 and conclusions end the paper. 

2. Stochastic Integrals 

Consider a family of increasingly fine triangulations schemes whereby the space dimension is 
discretized in multiples of an elementary space step hxm = Lx2~"^, m G N and we are interested 
in the limit as m ^ cxd. Let Xm = hxra"^ H X and consider the sequence of operators 

(2.1) £™ = ^A;" + M^)V™. 

defined on the 2™+^-dimensional space of all periodic functions fm ■ X„i M, where 

(2.2) V;V(2:) = + ^^") ~ - ^^'") . 
and 

(2.3) A™ /(.) = fi^+^^'^^) + n--h^rn)-2f{x) 

These definitions also apply to the boundary points by periodicity. We assume that m > mo 
where toq is the least integer such that 

(r,.-. a^(a:) ^ \ii{x)\ 

^ ' 2/i2 ^ 2h 

for all m > mo and all x G X^. 



Consider the kernel Um{x,x';t) of equation (2.6), i.e. the solution of the (forward) equation 



(2.5) —Um{x,x';t)=C'^*Um{x,x';t) 
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(2.7) Sx^x-x') 



where the operator £™* acts on the x' coordinate and the following initial time condition is 
satisfied: 

(2.6) Um{x,x';0) = h~^6x^{x - x'). 

where 

1 if X = x' mod 2Lx 
otherwise. 

If f{x, x') : Xm X Xjn M, consider the uniform norm 

(2.8) ||/m||m,oo= sup \f{x,x')\ 

and the graph norm 

(2.9) ||.flU,£=||.f||co+||/:r,/IU+||/:r/||oo. 

In (Albanese 2007c) we show that if the coefficients cr{xY and ijl{x) are uniformly continuous, 
then the sequence of kernels Um is Cauchy with respect to the graph norm. If in addition 
coefficients are Holder continuous and cr^ e C*^"'"", fi G C'^"'"" and 

(2.10) 7 = min{2,fc<, + a^,fc^ + a^} >0 
then 

(2.11) \\Um{t) - Um'{t)\\mX< Chl^ 

for all m' > m > nio- A similar bound also holds for the kernels obtained with a fully explicit 
Euler scheme 

(2.12) u'^{x,x';t) = h-^{l + dtmC^)^^] {x,x';t). 
where dtm is so small that 

(2.13) min 1 + 6tmC"'{x, x) > 0. 
In fact, we have that 

(2.14) \\Um{t)-u'^,{mm,C<chl^ 

for some constant c > 0. 

Let hyn be a monotonously dcscreasing sequence such as hyn as n ^ oo. Also let Lyn 
be a monotonously increasing sequence such that Lyn — > oo as n — > oo. Let 

(2.15) Yn = (llyn'^ n [—Lyn, Lyn]- 

We assume again periodic boundary conditions in the y direction and identify the two extreme 
points of Yn- Let's consider two-dimensional processes described by a sequence of Markov gen- 
erators of the form 

(2.16) C"^'^{x, y; x' , y') = C"'{x, x') + Q"'"(a^. V, v') 
where x,x' G X„ and y,y' € Yn- We assume that 

(2.17) ^ Q"^'^\x,y;x',y')=0 

y'eY„ 

for all values of x,y,x', so that the marginals with respect to the first process are the same as 
under the dynamics given by i.e. 

(2.18) exp(t£'"'")(x,2/;x',2/') = exp(f£")(x;x') 

for all triples x,y,x'. We are interested in the joint kernel 

(2.19) Un,n{x, y; x', y'; t) = ^ exp (a;, y; x' , y') 
and its convergence properties in the limit as m,n ^ oo. 
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The case of general multi-dimensional diffusion processes is discussed in (Albanese and Jones 
2007a). Results in that paper however do not apply to the class of processes in which we are 
interested here. In this article, we consider the case of path dependent processes which are not 
driven by their own diffusive dynamics and whose generator is not strongly elliptic, but rather 
hypoelliptic. Because of this reason, the joint kernel may be singular and a separate treatment 
is required. 



Definition 1. Consider a sequence of Markov generators of the form in {2.16) describing a 
bivariate process. One says that the second process is Abelian with respect to the first if the 
matrices Q^'^t of elements 

(2.20) Q-.'^(j;,y') = Q'"'"(^,y;a:',?/') 

are mutually commuting, i.e. if 

y^.^^j \yc x,x' ^ ^ x'\x"' \ ^x,x^x",x'" ^x',x'"^x,x' 

for all x,x' ,x" , x'" £ X^.. 

Definition 2. An Abelian process is called non-resonant if there is a sequence of non-singular 
transformations Vmnip, v) which diagonalizes simultaneously the kernels Q™'^i for all pairs x, x' € 
X^, i.e. if 

(2.22) V„™(p,y)Q™'"(x,y;x',y')V„^(p',2/') = A„„(x,x')5p,p'. 

for all x,x' G X„i. The index p ranges on a set denoted with Yn which has the same cardinality 
as Yn and is called the inverse lattice ofY^. 

In Sections [3] g] and [5] we consider in detail the prototypical example of Abelian process, 
namely stochastic integrals. In this case, a triangulation for the joint generator is given by 

/ cr(a;)2 ^(x) \ . , , u \x f ' \a{x)hx 



C"^-'\x,y;x',y')^ 



\ 2/?2 7h 
/ a{x)^ _ fi{x) 

Wx 

_ a{x) 



xm 
2 



2hr 



5x^ [x' -X - /i„„)(5y„ yy' - y- 

Sx„, {x' - X-\- ha:m)5Y^ [ u' ~ V + 



a{x)h^ 



5x^{x' - x)5YAy' ~y) + 6(x)4.'V"(y,2/') 



(2.23) 
where 

(2.24) 

and 

(2.25) 



V^5(2/) 



^Y-„ (y - y') 



gjy+hyn) -g{y-hyn) 

I 1 if y — y' mod 2Ly 
I otherwise. 



Diagonalisation is simply achieved by Fourier transforms of kernel 
(2.26) Tn{x,p;x',y) = e-'yPSxx, 

where 



(2.27) 



p e K„ = (hynZ) n 



hy 



and h,, 



Partial Fourier transforms in the p variable block-diagonalize the joint generator, reducing it to 
the form 



(2.28) 



{TnC^^'-T~'){x,p-x',p') = £™'"(x,x';p)V 
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where C"^'"{x,x';p) is a one-parameter family of matrices indexed by x,x' £ Xm- An explicit 
expression for the joint generator is given in equation (5.6) below. The kernel itself can also be 
expressed by means of a Fourier transform as follows: 
(2.29) 

Umn{x,p;x' ,p';t) EE (j^^^M„„ (i) J^„) (x, p; x' , p') = - — - — Syjp ~ p') exp (tC"'""{p)] {x,x'), 
We are interested in the convergence properties of the sequence 

(2.30) u,n{x,x';p,t) ^ lim -^Umn{x,p;x' ,p';t) 

in the limit as m ^ oo. To establish convergence, it is necessary to work in the representation of 
the Fourier transformed kernel. In fact, while the Fourier transformed kernel converges in a very 
strong graph-norm and the limit is entire analytic in p and has the same smoothness properties 
in X, x' as the diffusion kernel, the joint kernel expressed with respect to space coordinates can be 
quite singular even in simple cases. If for instance the functions a(x) and h(x) in (1.2) are zero, 
then the joint distribution \J (x, y; x' ,y' ; t) for the pair (xj, yt) is concentrated on the line yt — 0. 
If a partial Fourier transform with respect to the y variable is taken, the result is constant as 
a function of the dual variable p. This is a simple example of entire analytic dependency in p 
which however translates into a rough delta type singularity in the coordinate representation. 

Let g(x, x'; z) be a complex valued function defined for x, x' € x Xm and for z e C such 
that |z|< -fC where K > Q. Consider the vmiform norm 

(2.31) \\g\\m.,K,oo^ sup |5(x,x';z)|. 

\z\< K 

The graph norm of order (m, K, C) is defined as follows: 

(2-32) ||5l|m,if,£= llffm||m,K,oo + ||^r5l|m,-fS',oo + ||^™*ff||m,if,oo- 

Theorem 1. Let's assume the function b{x) is integrahle and let 
(2.33) B{x) = / b{x')dx' 



he a primitive. If the coefficients (t(x)^, m(2;) (^f^d the functions a(x) and B{x) are uniformly 
continuous, then the sequence of Fourier transformed kernels Umit) can he extended hy analyticity 
in p to an entire operator valued function and is Cauchy with respect to the ahove graph-norm 
for all K > 0. If in addition these coefficients are Holder continuous so that G Qk„,a„ ^ 
fi e C'=^""a', a e C'^-"", B e C''^^"^ and 

(2.34) 7 = min{2, fc„ + a^,k^ + a^^ ,ka + ks + as} > 
then, for all K > there is a constant c{K) such that 

(2.35) ||u™(i) - u„^'it)\\^^KX< c{K)hl^. 
for all m' > m > mg . 

Corollary 1. The limit kernel 

(2.36) u{x,x';z)= lim Um{x,x';z) 

is an entire analytic function of z which is in the domain of the operators and C* and all of 
their powers. The Fourier transformed kernel 

(2.37) u{x,x';y)= f e-*™w(x, x'; p) 



27r 

exists in the distribution sense for each fixed pair x, x' € [—Lx, Lx]. 
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Theorem 2. Let us assume that the coefficients a{x)'^ , and the functions a{x) and B{x) 

are uniformly continuous and fix a K > 0. Consider the kernels obtained with a fully explicit 
Euler scheme, i.e. 

, r f 1 

(2.38) u'J,{x,x';z,t) = — (l + <5t„£'"(^)) ^""^ {x,x'), 

where z g C : \z\< K. Assume that Stm satisfies the Courant condition 



(2.39) min ^[ 1 + StmC"^{x,x; z) ] > 
for all \z\< K. Then, there is a constant c{K) such that 

(2.40) \\uUt) - utmrn,K,C< c{K)hl^ 

for all m! > m> rriQ . 

This two Theorems are proved in the next two sections. An equivalent statement in a different 
representation which leads to fto's Lemma, the Cameron-Marin-Girsanov and the Feynman-Kac 
formulas is given in Section [Sj 

3. The Renormalization Group Argument 

The proof of convergence in graph-norm is based on a path-integral representation of the 
probability kernel. More precisely, let us defines a symbolic path 7 = {70, 7i, 72j .•.■} Si,s an 
infinite sequence of sites in Xm such that 7^ = 7j_i ± 1 for all j = 1, .... Let Fm be the set of all 
symbolic paths in X^,. Then the propagator admits the following representation: 

(3.1) Ujn{x,x';t) ^- — ^ ^ dsi ds2... dSqp{si,S2, ■..,Sq-i;j) 
where 

q-l 

(3.2) p(si,S2,...,s,„i;7)=e-^^^'"^^°'^''^n^™(Tj-i'Tj>~^''^'"'''^^"^^"^'^^«i-^«9- 

and Sq — t. 
Let 

7 € Tm and consider a path A'm(.; ti, ...tq\ 7) : — > X^n, left continuous and with right 
limits, taking the values {70,71,72,....} consecutively with jumps occurring at times tj,j = 
1, ..q — 1, with Q < ti < t2 < .. < ig-i < t. Let a{x) and h{x) be two uniformly continuous 
functions in AT = [—L^ , L^] and consider the integral 



I{ti,...tq_i;j) 



L 



o(A„(s;ti, ...^9:7)) ^^"*''^'^^^ ■■■^?'7) _|_ b{Xrn{s;ti, ...tq;j)) 



ds 



g-l q 

= X! «(7j)(7i+i - Ij) + X! ~ ^j)- 

(3.3) 

We are interested in evaluating the joint density 



Um{x,y;x' ,y';t) =- — ^ ^ dsi ds2... dsq 

p{si,S2, Sq_i;'y)S{l{si, ...Sq_i; 7) - (y' - y)) 



(3.4) 



8 



CLAUDIO ALBANESE 



We do so by expanding its Fourier transform 
(3.5) 

Urn{x,x';p;t) ^- — dsi ds2-. dsgp(si,S2, ...,s,-i;7)e'''^('*i- "«-i'^) 
The partial Fourier transform of the joint generator in (2.231 in the hmit as n — > oo is given 

by 

C"'{x,x';p) = hm 

^rl'-^ool^f (^,0;x,y,t)e .^Py=^— + —je ^™ ^ 



(3.6) 

where 



(3.7) cr{x,p, hxmY = cr(a;)^ cos(/ia;ma(a;)p) - isiii{hxmPa{x)) n{x) 

(3.8) /i(a;,P, /ixm) = cos{hxma{x)p) - ^(■^-)2^ 

(3.9) K(a:^,P, ft-xm) = _^ ^^"^ilT'xmPa{x)) ^^^^ _^ a{x)'^ {cos{h„na{x)p) - 1) - ipb{x). 

Notice that these three functions admit continuations as entire analytic functions in p. Let 
us denote them with a{x, z,hxm)^ , ^J'{x, z,hxm)j k{x, z,hxm)- Let L"^{x,x';z) be the analytic 
continuation of the operator above and consider the Fourier transformed joint propagator defined 
as follows: 



Um{x,x';z;t)^- — ^ ^ dsi ds2... dsqp{si, 

q=l^eT^:yn=x,^„=x' "1 "l-l 



S2,:;Sq-i;Z; 7) 



9=1 7er„:7o=a;,75 = 

(3.10) 
where 

p{si,S2....,s,.i-z-^)= — e-''''^^''°-'''-''^\{t^{^,.i,^f,z)^^^^^ 

^xm ■ -, 

(3.11) 

and Sq — t. 

If A' > 0, let us introduce the constants 

(3.12) So = inf a{x,0,h„n), M{K) = sup \fi{x, z,h„n)\ 

X e Am 

\z\< K 

and 

(3.13) "^liK) = sup ^|CT(a:, z, h^mW+Wix, z, hxm) - cr{x, 0, h^mW+h^nlKx, Z, Km)\- 

X e 

\z\< K 

Due to our assumptions, Sq > and T,i{K),M{K) < oo. 
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Figure 1. Contour of integration for the integral in (3.661. C+ is the countour 



joining the point D to the points E, A, B. C_ is the countour joining the point 
B to C to D. 



A symbolic path 7 — {70, 7i, 72, ■•■•} is an infinite sequence of sites in X^^ such that 7^ ^ Ij-i 
for all i — 1, .... Let r„ be the set of all symbolic paths in X^- The kernel of the diffusion 
process admits the following representation in terms of a summation over symbolic paths: 

C30 

(3.14) M„(x,y;z;t)=^2-« Yl W^{-/,q,z,t) 

''^^ J eT^ : Jo = x^-fq = y 
-Ij-ih 1 V j > 1 

where 
(3.15) 

"■xm Jo J si Js<j_i V 

with So = 0. 

Let us introduce the following Green's function: 

(3.16) Gm{x,y;z,u;) ^ Um.{x, y; z, t)e~"^*dt ^ h'^^ —{x,y)- 

Jo L"^(z) + iuj 

The Fourier transformed propagator can be re-obtained from the Green's function by evaluating 
the following contour integral: 

f djLO f did 

(3.17) u„,{x,y;z,t)= -^G,„(a;, y; z, w)e'"* + / ^G„,{x,y; z,uj)e"^* . 

Jc- Jc+ ^'^ 

Here, C+ is the contour joining the point D to the points E, A, B in Fig. [T| while C_ is the 
contour joining the point _B to C to D. The intent of the design of the contour of integration is 
to ensure that each point ui on the upper path C+ is separated from the spectrum of C 
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Lemma 1. For all K > 0, there are an integer m{K) and a constant c{K) such that 

duj 



(3.18) 



27r 



for all z such that \z\< K. 

Proof. The proof is based on the geometric series expansion 
1 °° \ 



C-^{z) + iu 



V2 / 5Cr(0, hxr.i)''^'' 



(3.19) 



1 .J 



Here (t(z, hxm)'^ and /i(z, /ia;m) are the multiphcation operators by cr(x, z, hxmY' and /i(a;, z, hxm), 
respectively. Convergence for w G C+ can be established by means of a Kato-Rellich type relative 
bound, see (Kato 1966). More precisely, for any a > 0, one can find a /3 > such that the 
operators V™ and A™ satisfy the following relative bound estimate: 

(3.20) ||V'"/||2<a||A"/||2+/3||/||2. 

for all periodic functions / and all m > mo. This bound can be derived by observing that V" 
and A™ can be diagonalized simultaneously by a Fourier transform and by observing that for 
any a > 0, one can find a /3 > such that 

1 



(3.21) 



< a 



for all m > mo and all k g Bm. 

Under the same conditions, we also have that 



(a(z, h^^Y - a(0, /i,™)^)A" + ^i{z, /i„„)V" / 



(3.22) 
Hence 



-a(0)^A"7 



lUJ 



< 



1^2^™^ 

2 ia^A 



1 



lU! 



— / 



< 1 



(3.23) 

where the last inequality holds if w € C+, if a is chos en su fficiently small and if m is large 
enough. In this case, the geometric series expansion in (3.19) converges in operator norm. 



The uniform norm of the kernel \Gm{x,y;uj)\ is pointwise bounded from above by h^^. 

Since the points B and D have imaginary part of height 4 ^ tJ^g integral over the contour 

C-|_ converges also and is bounded from above by cft,^ in uniform norm. 

□ 

Lemma 2. For all K > 0, if q > ^{2Y.i{KY + M{K)) we have that 



(3.24) 



W^™(7,(7;z,t) < 



■ exp 
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Proof. Let us define the function 

(3.25) 0(i) = ^^e l(t>0) 

where l(i > 0) is the characteristic function of M+. We have that 

(3.26) Wr„{i,q\z,t)<4>*'i{t) 

for all z such that |z|< K, where is the q— th convolution power, i.e. the g— fold convolution 
product of the function (j) by itself. The Fourier transform of </)(i) is given by 

The convolution power is given by the following inverse Fourier transform: 
(3.28) rm=l 0(-)V-^=(^) y_(l + ^) e-^^. 



^0 



Introducing the new variable z — 1 + ^^"2"' ; the integral can be recast as follows 

(3,2!.) r'it} ^ ""^Sr'" i„ '"•'■'^ (§'-'-'>) * 

where C^, is the contour in Fig. |2] Using the residue theorem and noticing that the only pole of 
the integrand is at z = 0, we find 

Making use of Stirling's formula q\ « \/2TTq'^^2e we find 



(3.31) (/.-^'(i) J^exp ( + glog + g(l - logg) 



If logq > log ^^2h^ * + ^' then we arrive at the bound in (|3.24|) 



□ 

Let us fix a if > 0. To prove the theorem, it suffices to consider the case m' — m + 1 for all 
values of m above mo. In fact, given this particular case, the general statement can be derived 
with an iterative argument. To this end, we introduce a renormalization group transformation 
based on the notion of decorating path. 

Definition 3. (Decorating Paths.) Let m > toq and let 7 = {yo, J/i, 2/2, .-.•} be o, symbolic 
sequence in Tm- A decorating path around 7 is defined as a symbolic sequence 7' = {j/Oi J/ij 2/2; ■■■■} 
with y[ e /i„j+iZ containing the sequence j as a subset and such that if y'j = yi and y'^, = J/i+i, 
then all elements y'^ with j < n < k are such that jy^ — yj\< hm+i- Let 'D„i+i{l) be the set of 
all decorating sequences around 7. The decorated weights are defined as follows: 



00 



(3.32) W^{^,q;z;t)^Yl E W^+^{i,q'-z-t). 

«-9 7' el?™+i(7) 

iq' = Iq 

Let us notice that these weights are not positive or even real values unless z — Q, but rather 
the depend analytically on the variable z. Finally, let us introduce also the following Fourier 
transform: 

poo ^ 1*00 

(3.33) W^{j,q;z;uj)= W^{j,q; z;t)e'^*dt, Wm{l,q; z;u;) = W„,{-/,q; z;t)e'^'dt. 
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Figure 2. Contour of integration Cb. for the integral in (3.29) 



Notation 1. In the following, we set h = h^.m+i so that hxm = 2ft,. We also use the Landau 
notation 0{h^) to indicate a Junction f{h) such that h^"f(h) is bounded in a neighborhood of 
0. 

Lemma 3. Let x,y E X„i and let C_ be an integration contour as in Fig. [7] Then 

, duj 



(3.34) 

Proof. We have that 
(3.35) 



2Gm+i(x, y] z; uj) - Gm{x, y; z; lo) e 



C- 



27r 



Oik'). 



2Gm+i{x, y; z; uj) - G„i{x, y; z;uj) ^^^2 « 

■'"^ 7 € r„ : 70 = a;,7g = y 
\lj -lj-i\= IVj > 1 

The number of paths over which the summation is extended is 

(3.36) N{^,q;x,y) = tt{7 g r,„ ■.j„ = x,jq = y, {jj - 7i-i|= IVj > 1} = (^l^k) 

where k — . Applying Stirling's formula we find 



^ (2Wr,r{-/,q;z;uj) -Wmh,q;z;uj) 



(3.37) 



^7 < 2" 
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Hence 



2Gm+i{x, y; z; - Gm{x, y; z; u) e 



< 



E 



max 



^ q=i V a g : 7o = 7, = y 
|7,--7,-i|= 1VJ>1 



5 du) 
2tt 



'2'Wmi'y, 9; z, uj) - Wm.i'y, q; z;uj)]e 



, duj 
2^ 



(3.38) 

for some constant c « > 0. It suffices to extend the summation over q only up to 



(3.39) 



2/i2 



To resum beyond this threshold, one can use the previous lemma. More precisely, we have that 

,„/ duj 



2Gm+i{x, y; z; oj) - Gm{x, y; z;uj)]e 



277 



< 



Cy/qn 



h 



'^'Wm{l,q\Z,Uj) - Wm{l,q\Z,UJl) 6 



27r 



max 

g, 7 G r„ : 7o = x, 7, = y 
|7,--7i-i|=lVi>l 

(3.40) 

Let us introduced the following abbreviated notations: 

(3.41) v{x, z) = a{x, z; ^h)"^ 

(3.42) m{x, z) = iJ,{x, z; ,h) 

(3.43) k(x, z) = k(x, z; ,h) 

where for the sake of keeping formulas short, we omit to denote the h = hx,m+i dependencies. 
To evaluate the resummed weight function, let us form the matrix 

/ v{x+h,z) j_ ^ j_ I. ^\ v{x+h,z) m{x+h,z) n. 

S2 ; -\- K{X-t n, Z) 2712 2h " 

v{x,z) . / \ v(x.z) m{x.z) 

P bK[X,Z) 2ft~ 

p v{x-h,z) m{x-h,z) v(x-h,z) / , x 



(3.44) jC{x,z) 
and decompose it as follows 



v{x,z) I m{x,z) 
2/i2 "I 2h~ 



\ 



(3.45) 
where 

(3.46) 



(3.47) 

(3.48) 
L-2.{x,z) 

and 
(3.49) 



C,{x, z) = -^£o{x, z) + j^'^iix, z) + C2{x, z) + hCi{x, z). 



^—v{x,z) ^v{x,z) 
jCo{x,z)= I ^v{x,z) —v{x,z) ^v{x,z) 







hv{x,z) —v{x,z). 



jCi{x,z) = \m{x,z) 



'-V--v{x,z) l\/fv{x,z)-lmix,z) 



— ^m{x, z) 



■\VZ^v{x,z) + \m{x,z) Vfv{x,z)^ 



\l^-^v{x,z) + k{x,z) iA™w(a;,0) - |V^m(a;,^) 

k{x, z) 

4 X "\-^^ 2 







^A'^'v{x,z) - ^V^m{x,z) 



■^A^v{x,z) + k{x,z), 



^W^k{x,z) -lAfm{x,z) 
C3{x,z) = I 

\A^mix,z) -V™fc(x,z), 
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Recall that all functions and operators above depend also on h = hx^m+i- 
Let us introduce the sign variable t — ±1, the functions 

(3.50) 0o(i, X, z, t) = 2C^{x, x + 2r/i; z)e*'^"(^'^'^)l(t > 0) 

(3.51) 0i(i,x,z,T) = 2£;"+i(a; + T/i,a; + 2T/i;z)e*^"(^'^)(x,a; + T/i)l(t > 0) 
and their Fourier transforms 

7 / X fv(x,z) m{x,z)\ /v{x,z) , . ^ 

Mu^, X, r) . A_i + )[-^- ^) + 

7/ /^ f(a:,z) m(a;, z) + V;"t;(a:, z) A^t;(a;, z) + V^m(x, z) Ag^m(x, z) , ,^.,2n 
0i(cj,a;,z,T) = ( — p hr \ 1 Th + U[h ) 



(3.52) 



where 



< x\ (—C{x, z) + iuj) ^ \x + rh > 



fO\ ( Sr.l 

(3.53) |a;>= 1 , and \x + rh >= \ 
We also require the functions 

(3.54) Mt, X, z) = et^"'('''^-'h{t > 0), Vi(^, = e'^'^y-'''> {x, x)l{t > 0) 
and the corresponding Fourier transforms 

^ / ^(^ 2^) \ ^ — ^ 

(3.55) ^Pq{uj , X , z) — I \-iu!j , ifji{uj,x, z) —< x\{—£{x, z) + ioj) \x > . 

If 7 is a symbolic sequence, then 

9-1 ^ 

(3.56) Wm{j,q;z,uj) = ■ipo{u;,'yq, z) J| 0o(w; ^, sgn(7j+i - 7^-)) 

j=o 
9-1 ^ 

(3.57) Wm(7,q;z;a;) = ■0i(w,75,z) J| <^i(a;; 7j,p,sgn(7j+i -7^))- 

j=o 

Let us estimate the difference between the functions (t>i{Lu, x, z, r) and (f>2{^, x, z, r) assuming 
that oj is in the contour C_ in Fig. [2j Retaining only terms up to order up to 0{h^), we find 

7 / N 1 2m(x,z)Th 4/i^ 

0o(w, a;, z, r) =1 H ^ r h ^ ^(k{x, z) - iuj)+ 

v{x,z) V(X,Z) 



4 



8m(x, z)Th^ . . . . 16h , , ,0 
\ \„ {k{x, z) - luj) + — r7,{k{x, z) - iujY + 0{h^). 

v{x,zy v{x,zy 

(3.58) 

A lengthy but straightforward calculation which is best carried out using a symbolic manipulation 
program, gives 

(w, x, z, r) = 1 + ; ; + r {k{x, z) - luj) - [8m{x, z) ~ V™z;(2;, z) 

v[x,z) v(x,z) ■'w(a;, z)^ 

+ R(x, z) ■ h^T + h'^Poix, z) + iLoK^PAx, z) + 0{h^) 

v[x, z)^ 



(3.59) 
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where 



R 



Po 



Pi 



(A™m + 2V;"/c)i;2 + (l6m - 2\'^v)kv 



{2m - 2V;"i;) V™fc + Uk' + (-4V^m - A^v)k ]v + [ - W + 2mV^i; + 2(V:,w)'' k 



28k - Am^ + 2mVZ'v - AvV^m - vA^v + 2(V™w)^ 



(3.60) 



For simplicity, we are not denoting here the dependency of all functions on (x, z). We have that 

( log0o(w;7i,^,sgn(7j+i - 7^-)) - log^i(w; 7^-, 2,sgn(7j+i -7^))) 
3=0 ^ ^ 



q~l 

E 



iwV;'u(7j,z) 



t;(7j,z)2 



+ R{-rj,z) )hhgn{jj + i - 7,) + {\\PQ\\K.oo + M\Pl\\K,oo+2\uj\^\\v-^\\K.,oo)0ih'^q) 

1 1 



^iioh^^^- ) +/i'(G(7„^)-G(7o,z)) + (||Polk,oo+kll|Pilk.oo+2|c^ni«-'|k,oo)0(/i49) 

(3.61) 



where G{x, z) is a primitive of R(x, z), i.e. 
(3.62) G(a;,z) = 



R{x',z)dx'. 



We conclude that there is a constant c(if ) > such that 
(3.63) 

J [H '^o('^' ^3'^' sgn(7j+i - 7j)) ~ Yi M^; sgn(7j+i - 7^)) j 



J=0 



2^ 



< c{K)K' 



for all g < (jmax- Here we use the decay of e*"* in the upper half of the complex lo plane to offset 
the u> dependencies in the integrand. Similar calculations lead to the following expansions: 
(3.64) 

iI}q{u),x,z) = — - + 0(a;/i^), i>i{uj,x,z) 



v{x, z) 

Since q < ch~^ and uj < |log/i|, we find 



, +0{iuh*)=lM^,x,z) + 0{iuh^). 
v(x, z) I 



(3.65) 



2G„+i(a:,y;z,[j) - G ,n{x , y; z , uj) e 



, doj 
2^ 



h 



By differentiating with respect to time in equation (3.661, we find that 



(3.66) ^u,n{x,y\t)= [ iwG„(a;,y;z;w)e*'^*^ + / ztjG,„(a;, y; z; w)e'"*^. 

dt Jc_ 2tx 7c+ 27r 

All the derivations above carry through. We conclude that for some constant c\{K) > and all 
\z\< K we have that 



(3.67) 



iujG,n{x,y]z;u))e" 
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and also 
(3.68) 



iuj[ 2Grn+iix,y;z;Lj) - Gm{x,y; z;u;) e 



27r 



< c 



Hence, the first time derivatives of the kernel satisfy a similar Cauchy convergence condition as 
the kernel itself. 

Let us notice that, if m is large enough, then the operator £™(z) is relatively bounded with 
respect to £™ in the following sense: 



(3.69) 



||/:/||<a(if)||£(z)/||+/3(i^)||/|| 



for some a{K), I3{K) > and all \z\< K. By inspection of the terms in C{z) one sees that 
such bound holds if h„n is small enough, i.e. for m larger than a threshold depending on K. K 
similar relative bound holds for the adjunct operator £* with respect to C{z)* . This concludes 
the proof of Theorem [l] 

□ 



4. Explicit Euler Scheme 

In this section we prove Theorem [2] A path- wise expansion for the time-discretization of the 
Fouier transformed kernel has the form 



9=1 7Sr„,:7o=a:,75=y fci = l k2=ki + l fcg=feg_i+l 

(4.1) {l + 6tt„,{^o, 7o; z)^ {6tY \{ £„(7,_i, 7, ; z) {l + ,5t£„(7, , 7,; z)^ 



where tq^i — t and fc^+i = N . In this case, the propagator can be expressed through a Fourier 
integral as follows: 

(4.2) u':^{x,y;z,t)= G^(x, y; z, c.)e-*^ 
where 

00 

(4.3) Gt{x,y-z,uj)=5tY,ut{x,y,z,j5t)e-''^^'K 

j=o 

The Fourier transformed propagator can also be represented as the limit 

St/„ „,. „ i\ _ 1™ f r^St/^ „,. „ , 



(4.4) <(x,?/;z,i) lim / GZ{x,y; z,uj)e 

where Ch is the contour in Fig. |3] This is due to the fact that the integral along the segments BC 
and DA are the negative of each other, while the integral over CD tends to zero exponentially 
fast as 3(w) — > cxD, where 9(a;) is the imaginary part of uj. Using Cauchy's theorem, the contour 
in Fig. |3]can be deformed into the contour in Fig. [l] To estimate the discrepancy between the 
time-discretized kernel and the continuous time one, one can thus compare the Green's function 
along such contour. Again, the only arc that requires detailed attention is the arc BCD, as the 
integral over rest of the contour of integration can be bounded from above as in the previous 
section. 

Let h = hxm and let us introduce the two functions 

(4.5) (l)o{t, X, z, r) EE 2£„(x, x + rh; z)e*^'"^"^^'''^h(t > 0), 

(4.6) (j)st{j,x,z,T) = 2Ljn{x,x + Th;z){l + 5tC„i{x,x;z)y \ 



iH-7r/6t 
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Figure 3. Contour of integration for the integral in (4.4) 



and the corresponding Fourier transforms 
(4.7) 

Xt f°° , f, -tujtd^^ fv{x,z,h) m{x,z,h)\ f v{x,z,h) 

<?>o(w,a:,z,T) = / (l)o[t,x,z,T)e 



27r 



/l2 



(4.8) 



w(a;, z 



We have that 



'v{x,z,h) m{x,z,h)\ ( . v{x,z,h) uj^ o- 
(/)5t(w, X, z,t)= [ — h T 1 H p -St + 0{St , 



(4.9) 



$o{u,x,z,t) 



2v{x, z, /i) 

where the last step uses the fact that St = 0{h?) 
Let us also introduce the functions 



h^5t + 0{h^6t^). = 0o(w, X, z, t) + 0{h'^), 



(4.10) Vo(t, a;, z, r) = e*^'"^^'^'^) l(i > 0), i^stU, r) = £ (l + Ji£™(x, x; z)) 



and the corresponding Fourier transforms 
(4.11) 

'v{x,z,h) 

tpoiuJ , X , Z , T) = I hlCJ 



k=l 



v{x,z,h)\ ^ 



/l2 



Again we find that 
(4.12) 



V'o(w, X, z, t) = -^if (w, x, z, t) + 0{h^). 
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If 7 is a symbolic sequence, then let us set 

9-1 ^ 

(4.13) W^h, q; z, uj) = '0o(w, 7,) Y[ <^o(w; ■^j,z, sgn(7j+i - 7^)) 

3=0 
«-i ^ 

(4.14) W^{-j,q; z,u) = -054(0;, 7,) J| 05t(w; 7^-, z, sgn(7j+i -7^))- 

j=o 

We have that 
(4.15) 

G™(x,y;z,w) - G™(a;,?/;z,w) =-^2^« ^ (m^^*(7, g; z, w) - #^(7, g; z, w 

■^"^ 7 e r„i : 70 = a;,7g ?/ 
l7j -7i-ih IVj > 1 

The integration over the contour in Fig. [T] can again be split into an integration over the 
countour C_ and an integration over C-|_. The integral over C+ can be bounded from above 
thanks to Lemma [T] Furthermore, we have that 

dw 



C- 



Gmi.^,y;z,uj) - Gm{x,y;z,uj) )e 



27r 



< 



Ch ^/On 



max 



g, 7 e r„ : 70 = x, 7, = y 
l7j ~7j-ih IVj > 1 



f duj 
2^ 



(4.16) 



To bound the time derivative, we have to consider 



o5t ^ 

'IT 



G'm(a;, y; z, w) - iLoGm{x, y; z, uj) ] e 



f duj 
2n 



(4.17) 



But, since St = 0{h^), also this difference is 0{h?). 



5. Convergence in the Ito Representation 
The convergence Theorem [T] admits a second formulation. Let us introduce the function 
(5.1) (l)7n{x) = ^ a{x'). 



We have that 



Jp<jim{x) 



K^'^xm '^il'xm/ 



( ^{xf _ fJ-jx) \ ^h^,„a{x)p. _ '^i^Y' 



Sx'a 



-(j{x,p,h^m) A"'{x,x') + fi{x,p, h„n)'^"'{x,x') 
- ^ '^™''^!^'"^°^^^V (a^) + ^(cos(/i„„a(a;)p) - 1) 

2 2 



(5.2) 



STOCHASTIC INTEGRALS AND ABELIAN PROCESSES 



19 



where a{x,p, hxm)'^ and fi{x,p,hxm) are defined in (3.7) and (3.8l and 



(5.3) 

Hence, we have that 



Vj/(x) 



f{x+h„a) - f{x) 



(5.4) 
where 
(5.5) 
Here 



a{xf 



U{x,p, h^m) = ipb{x) + ^ (^/i-^V+e-^''^'"'^(-)P^ , 



r^{x,x';p)=(^h^^Wte~^'^-<^'>P^ ( 
(5.6) 



'"^ + ^^)v^(.,.') + ^4 



Jp{(t>m{x')--ip't>m(x)) 



and 
(5.7) 



V™+/(a:) = 



f{x + Km) - fjx) 



The operator U^{x,x';p) is equivalent to C"^{x,x';p) up to a non-singular linear transfor- 
mation. We say that L"^(x,x';p) is the Fourier transformed generator in the Ito representation. 
Also notice that the weak limit of this operator as m — s- cx) is equal to 

,5.8) = ^£ - 

The Fourier transformed kernel is given by 

Urn{x,x';p) = e^f(*"(-)-^'"(-'))[/„(x,a;';p) 

(5.9) 
where 

Um{x,x';p) = /i^^e*-^'"(a;,a:')- 

(5.10) 

The joint kernel is thus given by the following formula whose continuum analog was found in 
(Girsanov 1960), (Cameron and Martin 1949), (Feynman 1948) and (Ito 1949): 

dp 



u{x, I;x ,1 ;t) 



2tt 



.jpli'-i-^^A^')+M^)]U^^x,x'-p). 



(5.11) 

Theorem 3. Under the same assumptions of Theorem^ for all K > there is a constant c{K) 
such that 

(5.12) \\Uiit)-UUt)\L,K,c<^iK)hZm. 

for all m' > m > uiq. A similar bound also holds for the kernels obtained with a fully explicit 
Euler scheme, i.e. 

(5.13) U^ix,x';p,t) = ^(l-K5t™L™(p))[^] ix,x'), 
where 5t„i is so small that 

(5.14) min I + 5tmt^{x,x;p) > Q. 
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In this case, there is a constant c{K) such that 

\m,K,C 



(5.15) \\Um{t) - U^mm.K.C< ciK)hl, 



for all m' > m > mg . 

It is possible to retrace the argument in the previous section except for replacing 

(5.16) K.{x,p, h^rn) Cm{x,p, h.^yn) + hx^rm{x,p, kxm) 



(6.3) A"\x,y;x',y') 



where k{x,p, h^m) is the function in (3.9 1. All arguments go through unaltered under the same 
conditions. As a consequence we conclude that the kernels 

(5.17) exp ( + Cm{z) + Kn.rm{z))){x, x') 

converge in graph norm, uniformly on discs z d C : \z\< K , for any K > 

6. Other Abelian Processes 

In this Section, we give to examples of Abelian processes which emerge from applications and 
are not stochastic integrals. 

6.1. The Sup Process. Consider the sup process 

(6.1) yt = sup Xg. 

e[o,t] 

the sup is always attained as an element of the underlying space Xm, it is natural in this case to 
restrict the attention to the case where n — m and = X„i. The joint generator is given by 

y- x', y') x')A^{x, y; x', y') 

(6.2) 

where x, y, a;', y' € X^a and we set 

5yyi if x' < y 
5x'y' if x' > y. 
Consider the matrix 

(6.4) VUx, y; x\ y') = 5{x - x')l{y' > y) 
and its inverse 

(6.5) V;;^(a;,y;a;',y') = 6{x - x){S{y' - y) - S{y' - y~ hr^^)). 
Consider the one parameter family of operators C"^{y) such that 

(6.6) {V-,'£^V„^){x,y;x',y')^£'^{x,x';y)Syy,. 
We have that 

(6.7) C^{x,x';y) = C^{x;x')l{x' <y). 

Hence, C"^(x,x';y) is the Markov generator of the underlying process with absorption in the 
interval [y, L^,]. 

Let us notice that the kernel can be obtained as follows: 

(6.8) ■Um(t) = Vm{tm(t)V,7/ where Urn{x,y-x' ,y' \t) = 5{y - y')eiiY> {tC"\y)^ {x,x'). 

A more explicit way of expressing the joint kernel is 
(6.9) 

.„(., t) = 5iy - ,')^™((., V, t) + Ky' > .) (^™(x, t) - u^x, x'; y' - h.^, t)'^ 

Convergence in this case can be established along the same lines as done for stochastic in- 
tegrals. The situation is simpler in that only the consideration of the kernel itself, i.e. the 
p = z = case with the notations in Sections 2 and 3, is needed. The additional complication is 
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that we need to consider absorbing boundary conditions. This implies a few marginal changes 
to the derivation above, as when a path arrives to an absorbing lattice point, it stays constant 
thereafter. Since this applies both to paths and decorating paths in a finer lattice, upon arriving 
to an absorption point the dynamics is trivial in either case and the final bounds given still hold. 

6.2. Discrete Time Processes. This section is based on work in collaboration with Manlio 
Trovato (Albanese and Trovato 2005) and Paul Jones, see (Albanese and Jones 2007&). 

An important class of path-dependent options requires computing the joint distribution of 
the underlying lattice process and of a discrete sum of the following form; 

N 

(6.10) yt = ^ip{xt,_,,xti,ti) 

i=l 

where N is an integer, ti = iAT and T = NAT. Consider the elementary propagator 

(6.11) C/™(.xi,.X2) = r^e(™(xi,,T2). 

To find the joint transition probability, one can again discretize the variable yt in the lattice 
Yn = hyn'Z n [—Ly,Ly]. As opposcd to lifting the generator as done above for the other cases, 
here we lift the elementary propagator itself and form the joint propagator 

(6.12) Umn{xi,yi;X2,y2) = Um{xi,X2)6{y-i - J/2 + [lp{xi,X2)h~^]). 

This lifted operator can be block-diagonalized by means of a partial Fourier transform. Consider 
the Fourier transform operator (p) of matrix elements 

(6.13) C/„(xi,X2;p) = lim Umn{xi,yr,X2,y2)e-"'^y'-y'^ = Um{xi,X2)e-'^^^'''''''\ 

n^oo 

Then we have that 

(6.14) Jim {U^^){xuyi;x2,y2) = ^e'P^y^-y^\U^){xi,X2;p) 

Convergence in the graph-norm in this case descends directly from the convergence of the one- 
period kernel Um{xi,X2)- 

7. Conclusions 

We obtained bounds on convergence rates for explicit discretization schemes to the Fourier 
transform of joint kernels of one-dimensional diffusion equations with continuous coefficients and 
a stochastic integral. We consider both semi-discrete triangulations with continuous time and 
explicit Euler schemes with time step small enough for the method to be stable. The proof 
is constructive and based on a new technique of path conditioning for Markov chains and a 
renormalization group argument. Convergence rates depend on the degree of smoothness and 
Holder differentiability of the coefficients. The method also applies to a more general class of 
path dependent processes we call Abelian. Examples of Abelian processes beside stochastic 
integrals are the sup process and discrete time summations. 
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